Towards a new partially integrated transport model for coarse grid and unsteady turbulent flow simulations

A new two-equation model is proposed for large eddy simulations (LESs) using coarse grids. The modeled transport equations are obtained from a direct transposition of well-known statistical models by using multiscale spectrum splitting given by the filtering operation applied to the Navier–Stokes equations. The model formulation is compatible with the two extreme limits that are on one hand a direct numerical simulation and on the other hand a full statistical modeling. The characteristic length scale of subgrid turbulence is no longer given by the spatial discretization step size, but by the use of a dissipation equation. The proposed method is applied to a transposition of the well-known k- statistical model, but the same method can be developed for more advanced closures. This approach is intended to contribute to non-zonal hybrid models that bridge Reynolds-averaged Navier–Stokes (RANS) and LES, by using a continuous change rather than matching zones. The main novelty in the model is the derivation of a new equation for LES that is formally consistent with RANS when the filter width is very large. This approach is dedicated to applications to non-equilibrium turbulence and coarse grid simulations. An illustration is made of large eddy simulations of turbulence submitted to periodic forcing. The model is also an alternative approach to hybrid models. PACS 47.27.Eq     

后台阶流动的数值模拟

访述了大涡模拟的基本思想，指出大涡模拟的效率主要取决于四个因素，即流动中须有大尺度涡存在、合理的计算格式、合适的滤波器和亚格子应力模型。在深入考虑粘性不可压缩流Navier—Stokes方程各个子项作用的基础上，提出二阶全展开Euler—Taylor—Galerkin有限元方法作为大涡模拟的离散格式，并采用Gauss滤波器，对典型算例——后台阶处的流动进行大涡模拟，计算结果与相关文献符合的很好。从计算结果还可以看出大涡模拟与二阶全展开ETG有限元方法的结合在捕捉涡系及反映涡动时变过程方面具有明显的优势，说明大涡模拟适合于边界几何形状复杂区域流动的模拟。同时应用二阶全展开ETG有限元方法对低雷诺数粘性不可压缩后台阶流动进行了计算，得到与相关文献符合良好的计算结果，即该方法也可独立用于对低雷诺数粘性不可压缩流动的计算。     

Influence of Molecular Characteristics of a Substance on Detonation Phenomena

In the present study, the molecular-dynamics method has been used to examine the influence of the thermal effect of the chemical reaction on processes in a detonating molecular crystal (propagation velocity and structure of the detonation wave, delay time of the chemical reaction, and energy macrocharacteristics of the crystal in the disturbed region). Molecular-dynamics data are compared with predictions of the continuum theory of detonation; in particular, fulfillment of the Chapman–Jouguet condition is verified.     

亚临界雷诺数下圆柱绕流的大涡模拟

本文应用Smagorinsky涡粘性模式和二阶精度的有限体积法对圆柱绕流的流场进行大涡模拟．求解了非正交曲线坐标系下的N-S方程,对雷诺数为100和20000的工况进行了计算．计算结果与实验及动力涡粘性模式的结果进行了比较,表明计算对于层流及高亚临界雷诺数的湍流流动是合理的     

白洋淀内水体的流动特征研究

研究了风对白洋淀内水体循环流动的影响，研究指出：(1)风应力可看作是水体循环 的主要驱动力. (2)在风应力的作用下，水体表层的流动主方向与风的方向一致，下层水体 的流动方向与表层方向相逆或近似相逆. (3)表层以下水体流动复杂，次表层水平流动有明 显的双涡结构. (4)淀底地形对水体垂向流动影响很大，垂向速度变化复杂，垂向上有双涡 流动现象. 从动力学角度研究了淀内水体流动的情况.     

Coherent fine scale eddies in turbulence transition of spatially-developing mixing layer

To investigate the relationship between characteristics of the coherent fine scale eddy and a laminar–turbulent transition, a direct numerical simulation (DNS) of a spatially-developing turbulent mixing layer with Re_ω,0 = 700 was conducted. On the onset of the transition, strong coherent fine scale eddies appears in the mixing layer. The most expected value of maximum azimuthal velocity of the eddy is 2.0 times Kolmogorov velocity (u_k), and decreases to 1.2u_k, which is an asymptotic value in the fully-developed state, through the transition. The energy dissipation rate around the eddy is twice as high compared with that in the fully-developed state. However, the most expected diameter and eigenvalues ratio of strain rate acting on the coherent fine scale eddy are maintained to be 8 times Kolmogorov length (η) and :β:γ = −5:1:4 in the transition process. In addition to Kelvin–Helmholtz rollers, rib structures do not disappear in the transition process and are composed of lots of coherent fine scale eddies in the fully-developed state instead of a single eddy observed in early stage of the transition or in laminar flow.     

Near-Wall Modification of an Explicit Algebraic Reynolds Stress Model Using Elliptic Blending

The elliptic blending approach is used in order to modify an Explicit Algebraic Reynolds Stress Model so as to reproduce the correct near wall behaviour of the turbulent stresses. The anisotropy stress tensor is expressed as a linear combination of tensor bases whose coefficients are sensitised to the non-local wall-blocking effect through the elliptic blending parameter γ. This parameter is obtained from a separate elliptic equation. The model does not use the distance from the wall thus it can be easily applied to complex geometries. It is validated against detailed DNS data for mean and turbulence quantities for the case of flow and heat transfer between parallel flat plates at three Reynolds numbers as well as against experimental data for the flow in a backward facing step at Re _H = 28,000. The comparison with DNS results or experiments is quite satisfactory and shows the validity of the approach.     

Analysis of Unsteady Motion with Respect to Noise Sources in a Gas Turbine Combustor: Isothermal Flow Case

In order to evaluate the direct and indirect contributions to the total combustion noise emission, a combustion chamber consisting of a swirl burner and an exit nozzle of Laval-shape, representing a gas turbine combustor, is investigated by means of experiments and large eddy simulation. Focused on the isothermal flow case first and encouraged by a good overall agreement between the LES and the experimental data for the flow field, a first characterisation of the flow with respect to noise sources is performed. To analyse acoustic properties of the flow, time and length scales are evaluated inside the combustor. Furthermore, the evidence for the existence of a precessing vortex core (PVC), typical for configurations with swirl, is revealed. Finally, the effect of the PVC on the flow inside the Laval nozzle is discussed.     

Subgrid-modelling in LES of compressible flow

Subgrid-models for Large Eddy Simulation (LES) of compressible turbulent flow are tested for the three-dimensional mixing layer. For the turbulent stress tensor the recently developed dynamic mixed model yields reasonable results.A priori estimates of the subgrid-terms in the filtered energy equation show that the usually neglected pressure-dilatation and turbulent dissipation rate are as large as the commonly retained pressure-velocity subgrid-term. Models for all these terms are proposed: a similarity model for the pressure-dilatation, similarity andk-dependent models for the turbulent dissipation rate and a dynamic mixed model for the pressure-velocity subgrid-term. Actual LES demonstrates that for a low Mach number all subgrid-terms in the energy equation can be neglected, while for a moderate Mach number the effect of the modelled turbulent dissipation rate is larger than the combined effect of the other modelled subgrid-terms in the filtered energy equation.     

Optical Analysis and Qualitative Discrimination of Vortex-Flame Interaction in a Plane Premixed Shear Layer

To obtain practical schemes of vortex–flame interactions, a series of organized eddies formed in the plane premixed shear layer is investigated, instead of a single vortex ring or a single vortex tube. The plane premixed shear layer is first formed between two parallel uniform propane–air mixture streams. For getting clear qualitative pictures of vortex–flame interactions in the plane premixed shear layer, two extreme ignition points are assigned; one is assigned at the center of an organized eddy where the vortex motion plays an important role, the other at the midpoint between two adjacent organized eddies where the rolling-up motion prevails. A premixed flame is initiated by an electric discharge at one of the two assigned points and propagates either in the large scale organized eddy or along the interface between two uniform mixture streams. Propagation and deformation processes of the flame are observed using the simultaneously two-directional and high-speed Schlieren photography. The tangential velocity of organized eddy and the equivalence ratio of premixed shear flow are varied as two main parameters. The outline of propagating flame after the midpoint ignition is numerically analyzed by superposing the flame propagation having a constant burning velocity on the vortex flow field simulated with the discrete vortex method. The results obtained show that there exists another type of vortex–flame interaction in the plane shear layer in addition to the vortex bursting, and that it is caused by the rolling-up motion particular to the coherent structure in the plane shear layer and is simply named the vortex boosting. It is qualitatively concluded therefore that, in the ordinary turbulent premixed flames formed in the plane premixed shear layer, these two fundamental vortex-flame interactions get tangled with each other to augment the propagation velocity. An empirical expression which qualitatively takes into account of the effects of both vortex and chemical properties is finally proposed.